Dynamics of strangeness production in heavy-ion collisions near threshold energies 



O 

(N 
-i— > 
O 

O 

m 



o 

=5 



> 

o 

d 
o 



Zhao-Qing Fen^] and Gen-Ming Jin 
Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, People's Republic of China 

(Dated: October 14, 2010) 

Within the framework of the improved isospin dependent quantum molecular dynamics (ImlQMD) 
model, the dynamics of strangeness (K 0,+ , A and E _,0 ' + ) production in heavy-ion collisions near 
threshold energies is investigated systematically, in which the strange particles are considered to 
be mainly produced by the inelastic collisions of baryon-baryon and pion-baryon. The collisions in 
the region of supra-saturation densities of the dense baryonic matter formed in heavy-ion collisions 
dominate the yields of strangeness production. The total multiplicities as functions of incident 
energies and collision centralities are calculated with the Skyrme parameter SLy6. The excitation 
function of strangeness production is analyzed and also compared with the KaoS data for the K + 
production in the reactions 12 C+ 12 C and 197 Au+ 197 Au. 
PACS number (s): 25.75.-q, 13.75. Jz, 25.80. Dw 



Particle production in relativistic heavy-ion collisions 
has been investigated as a useful tool to extract the infor- 
mation of the nuclear equation of state (EoS) of isospin 
asymmetric nuclear matter under extreme conditions in 
terrestrial laboratories, such as high density, high tem- 
perature and large isospin asymmetry etc. Besides nu- 
cleonic observables such as rapidity distribution and flow 
of free nucleons or light clusters (such as deuteron, triton 
and alpha etc.), also mesons emitted from the reaction 
zone can be probes of the hot and dense nuclear matter. 
In particular, the strangeness productions (K° :+ , A and 
yj-,o,+-j m nea vy-ion collisions in the region of 1A GeV 
energies have been extensively investigated both exper- 
imentally P, [H and theoretically ■ Kaons (K° and 
K + ) as a probe of EoS are produced in the high den- 
sity phase without subsequent reabsorption process. It 
was noticed that the kaon yields are sensitive to the EoS 
in theoretical investigations by transport models. The 
available data already favored a soft EoS at high densi- 
ties. The K° /K + ratio was also proposed as a sensitive 
probe to constrain the high-density behavior of the sym- 
metry energy 0, @|. Up to now, the difference of the 
symmetry energy at supra-saturation densities predicted 
by transport models is huge. It is not only in understand- 
ing the reaction dynamics, the high-density behavior of 
the symmetry energy also has an important application 
in astrophysics, such as the structure of neutron star, the 
cooling of protoneutron stars, the nucleosynthesis during 
supernova explosion of massive stars etc [9]. With the 
establishment of new-generation radioactive beam facil- 
ities in the world, such as the CSR (IMP in Lanzhou, 
China), FAIR (GSI in Darmstadt, Germany), PJKEN 
(Japan), SPIRAL2 (GANIL in Caen, France) and FRIB 
(MSU, USA) , the high-density behavior of the symmetry 
energy is being studied more detail experimentally. Pre- 
cise description of the experimental data is still necessary 
by improving transport models or developing some new 
approaches, especially for constraining the high-density 
behavior of the nuclear symmetry energy. 

The ImlQMD model has been successfully applied 



to treat dynamics in heavy-ion fusion reactions near 
Coulomb barrier and also to describe the capture of two 
heavy colliding nuclides [H| • The ground state prop- 
erties of a single nuclide is described well in the ImlQMD 
model. Further improvements of the ImlQMD model 
have been performed in order to investigate the pion 
dynamics in heavy-ion collisions and to extract the in- 
formation of isospin asymmetric EoS at supra-saturation 
densities fl3 - ll4) . The strangeness productions at near 
threshold energies are to be included in inelastic baryon- 
baryon and pion-baryon collisions and investigated in this 
work. 

In the ImlQMD model, the time evolutions of the 
baryons and mesons in the system under the self- 
consistently generated mean-field are governed by Hamil- 
ton's equations of motion, which read as 



dH 



dH 
dpi' 



(1) 



Here we omit the shell correction part in the Hamiltonian 
H as described in Ref . . The Hamiltonian of baryons 
consists of the relativistic energy, the effective interaction 
potential energy and the momentum dependent part as 
follows: 



u„ 



(2) 



Here the p i and rrii represent the momentum and the 
mass of the baryons. The momentum dependent term is 
taken as the same form in Ref. }15| . which reduces the 
effective mass in nuclear medium. 

The effective interaction potential is composed of the 
Coulomb interaction and the local interaction 

Uint = Ucoul + Ul oc . (3) 

The Coulomb interaction potential is calculated by 

Ucoul = JE ^er/(r«/>/4L) (4) 

*>...,, Tij 
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where the ej is the charged number including protons 
and charged resonances. The = |r^ — Yj | is the relative 
distance of two charged particles. 

The local interaction potential energy is derived di- 
rectly from the Skyrme energy-density functional and ex- 
pressed as 



V loc {p{v))dr. 



(5) 



The local potential energy-density functional reads 



Vloc(P) =~ — 
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where the p n , p p and p = p n + p p are the neutron, 
proton and total densities, respectively, and the 8 — 
(p n — p p )/(p n + Pp) is the isospin asymmetry. Here, all 
the terms in the Skyrme energy functional are included 
in the model besides the spin-orbit coupling. The coef- 
ficients a, /3, 7, g S ur, glmn 9t are related to the Skyrme 
parameters to , ii , i 2 ,^3 and xq, x\, cc 2 , £3 as [1 li - 
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(7) 
(8) 
(9) 
(10) 



The parameters of the potential part in the bulk sym- 
metry energy term are also derived directly from Skyrme 
energy- density parameters as 
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Csym^-^^ 2 ) P0 /3 [3il^l-t2(5x 2 +4)]. (11) 

In this work, we use the Skyrme force Sly6 in the Im- 
IQMD calculations, which gives the modulus of incom- 
pressibility of symmetric nuclear matter about 230 MeV. 

Analogously to baryons, the evolution of mesons (here 
mainly pions and kaons) is also determined by the Hamil- 
tonian, which is given by 



N M 

H M =J2{ v P 0Ul +^p>)) 



(12) 



Here the Coulomb interaction is given by 
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where the Nm and Nb are the total numbers of mesons 
and baryons including charged resonances. Here, we use 
the energy in vacuum for mesons, namely cj(pj,pi) = 



vw 



with the momentum p, and the mass tum of 



the mesons. 

A hard core scattering in two particle collisions is as- 
sumed in the simulation of the collision processes by 
Monte Carlo procedures, in which the scattering of two 
particles is determined by a geometrical minimum dis- 
tance criterion d < ^/0.1a tot /7r fm weighted by the Pauli 
blocking of the final states [lg, [l7 1 . Here, the total cross 
section a t ot in mb is the sum of the elastic and all inelas- 
tic cross section. The probability reaching a channel in a 
collision is calculated by its contribution of the channel 
cross section to the total cross section as P c h = <J c h / &tot ■ 
The choice of the channel is done randomly by the weight 
of the probability. The primary products in nucleon- 
nucleon (NN) collisions at the 1A GeV energies are the 
resonances A(1232), N*(1440) and the pions. The reac- 
tion channels are given as follows: 

AAoAA, NN -H- AA*, NN <h> AA, 

A Att, N* O Att, AA -> AAtt(s - state). (14) 

The cross sections of each channel to produce resonances 
are parameterized by fitting the data calculated with the 
one-boson exchange model (H)]. In the 1 A GeV region, 
there are mostly A resonances which disintegrate into a 
7r and a nucleon, however, the A* yet gives considerable 
contribution to the high energetic pion yield. The en- 
ergy and momentum dependent decay width is used in 
the calculation [IJJ] . The strangeness is created by inelas- 
tic hadron-hadron collisions. We include the channels as 
follows: 

BB -> BYK, BB ->■ BBKK, Btt -> YK, Bit AAA, 
Ytt -» AA, AA ->■ Yn, YN -> AAA. (15) 

Here the B strands for (N, A) and Y(A, £), K(K°, K+) 
and K(K°, K~). The parameterized cross sections of 
each isospin channels BB — s- BYK [Tj[ are used in the 
calculation. We take the parametrizations of the chan- 
nels Btt -> YK besides the Nir -> AA reaction 
21| . The results are close to the experimental data at 
near threshold energies. The cross section of antikaon 
production in inelastic hadron-hadron collisions is taken 
as the same form of the parametrization used in hadron 
string dynamics (HSD) calculations (13]. Furthermore, 
the elastic channels are considered through the channels 
KB — > KB and KB —t KB and we use the parametriza- 
tions in Ref. 23 1. 
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FIG. 1. Time evolution of the central density, pion and strangeness produced in the reaction 197 Au+ 197 Au for head-on collisions 
at incident energy 1.5A GeV. 



We calculated the production of pions, kaons and hy- 
perons and compared with the evolutions of the central 
density in central 197 Au+ 197 Au collisions at incident en- 
ergy 1.5A GeV as shown in Fig. 1. One can see that 
the pions and strange particles are mainly produced at 
supra-saturation densities (p > po). The observables can 
be probes to extract the high-density information of EoS. 
The pion yields saturate at time of the order of 30 fm/c, 
and the kaons are produced at early stage of the reaction 
around 25 fm/c. The results are similar to the reports 
in Ref. [7| but more pronounced in the RBUU calcu- 
lations. The Bit — » YK contributes to be about 1/3 
of the total kaon production, which retards the profile 
of the kaon yields in the calculations. The saturation 
of the hyperons is more slowly than the cases of pions 
and kaons owing to the exchange reactions Yir — > NK 
and NK — > Y~k . The isospin difference is pronounced in 
the temporal evolutions, which enables the observables 
as probes of the high-density behavior of the symmetry 
energy, such as the ratios tt~ /ir + , K°/K + and E _ /E + . 

The strange particles in heavy-ion collisions are mainly 
created at supra-saturation densities. To extract the 
basic information of the compressed nuclear matter at 
high densities formed in relativistic heavy-ion collisions, 
the centrality and incident energy dependence of the 
strangeness production was calculated as shown in Fig. 
2 and in Fig. 3 for the 197 Au+ 197 Au reaction. A larger 
domain of the high-density phase diagram is formed in 
central collisions and with increasing incident energy, 
which enhances the probabilities of NN inelastic collisions 
and hence the strangeness production. The neutron- 
neutron collisions contribute the productions of K° and 
Therefore, the neutron-rich systems augment the 
values of the K°/K+ and IT /£+ ratios. 

The reliability of the calculations on the strangeness 
production can be checked from the available experimen- 



tal data. Shown in Fig. 4 is the calculated excitation 
functions of strange particles for the heavy 197 Au+ 197 
and the light 12 C+ 12 C reactions and compared with the 
KaoS data for the K + production [l(. The experimental 
data can be well reproduced by calculations besides at 
very low threshold energies owing to the limited statis- 
tics. The larger cross sections of the strangeness produc- 
tion for heavy systems result from the larger region of the 
high-density phase diagram. The difference of production 
of the same isospin particles such as K 0,+ or is de- 

viated from the isospin effects of collision systems. Influ- 
ence of the symmetry energy and the in-medium poten- 
tials on the strangeness production is being investigated. 

In summary, the production of strangeness in heavy- 
ion collisions in the region of 1A GeV energies is investi- 
gated within the framework of the ImlQMD model. The 
strange particles are produced at the supra-saturation 
densities and emitted at the early stage of the phase di- 
agram comparing with the pion production. Strangeness 
exchange process retards the saturation of A and £ pro- 
duction. Higher incident energy and central collisions 
enhance the domain of the high-density phase diagram, 
hence lead to the increase of the strangeness multiplicity. 
A larger high-density region of compressed nuclear mat- 
ter is formed in heavy collision systems comparing with 
the light systems, which increases the strangeness pro- 
duction. The experimental data of the K + production 
cross sections at near threshold energies can be repro- 
duced rather well. 
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toral Project of Chinese Academy of Sciences; and the 
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FIG. 3. The same as in Fig. 2, but the calculated excitation functions of strangeness production for head-on collisions. 




FIG. 4. Comparison of the calculated strangeness production cross section per mass number (er/A) and the KaoS data for the 
K+ production for the system 197 Au+ 197 Au and the reaction 12 C+ 12 C. 
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